Hot electron relaxation in a heavy fermion system with tuned doping. 
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Femtosecond time-resolved optical spectroscopy was used to systematically study photoexcited 
carrier relaxation dynamics in the intermediate-valence heavy fermion system Ybi-^Lu^Ala (0 < 
x < 1). Given the demonstrated sensitivity of this experimental technique to the presence of the 
low energy gaps in the charge excitation spectrum, the aim of this work was to study the effect 
of dilution of the Kondo lattice on its low energy electronic structure. The results imply that in 
Ybi-^Lu^Als the hybridization gap, resulting from hybridization of local moments and conduction 
electrons, persists up to 30% doping. Interestingly, below some characteristic, doping dependent 
temperature T*(x) the relaxation time divergence, governed by the relaxation bottleneck due to 
the presence of the indirect hybridization gap, is truncated. This observation is attributed to the 
competing ballistic transport of hot electrons out of the probed volume at low temperatures. The 
derived theoretical model accounts for both the functional form of relaxation dynamics below T*(x), 
as well as the doping dependence of the low temperature relaxation rate in Ybi-aUuaAbj. 



I. INTRODUCTION 

A reasonable understanding of the Kondo problem was 
achieved in the case of isolated Kondo impurities by the 
Anderson impurity model (AIM). Despite the fact that 
many heavy electron systems are stoichiometric with a 
local moment in each unit cell, many of their physical 
properties can be well accounted for by the AIM. On the 
other hand, some low-temperature properties such as op- 
tical conductivity are fundamentally different than those 
expected of the AIM. These are consistent with the renor- 
malized band behavior, captured by the Anderson lattice 
model. Here hybridization of local moments and conduc- 
tion electrons leads to the presence of the hybridization 
gap (HG) in the vicinity of the Fermi levelj 1 "^ where in 
heavy fermions the Fermi level lies close to the edge of 
the lower hybridized band.— ^ 

Recent experiments on the photoexcited carrier relax- 
ation dynamics in several heavy electron compounds have 
shown that the relaxation of the electronic system back 
to equilibrium is extremely sensitive to the underlying 
low energy electronic structure and strongly depends on 
the temperature (T) and excitation level^ 5 -^^ It has 
been shown, that the bottleneck in carrier relaxation 
is governed by the presence of the HG near the Fermi 
leveli^£ Both T and excitation density dependence of 
the relaxation process can be well accounted for by the 
phenomenological Rothwarf- Taylor (RT) modet 5 ^ ' 10 ^ 1 
which was originally developed to describe the relaxation 
in superconductors driven out of equilibrium ^ 11 ^ 11 We 
should note that the alterantive model presented in Refs. 
[6,7] is esentially also a bottleneck model. Here the en- 
ergy and momentum conservation law leads to the sup- 
pression of e-ph scattering near in the case of sound 
velocity exceeding the Fermi velocity. This suppression 
of e-ph scattering near E / effectivly acts as a gap in the 



DOS. However, since the temperature and excitation de- 
pendence of the relaxation dynamics in a Kondo insu- 
lator SmBg is nearly identical to that of heavy fermion 
systemsj 4 ^ 5 - it was naturally to assume that the origin of 
the bottleneck in the entire class of materials was in the 
presence of the indirect hybridization gap^i 5 - 

In this paper we report first systematic studies of pho- 
toexcited carrier relaxation dynamics in an intermediate- 
valence system Ybi-xLu^Ala. In Ybi-zLuzAla concen- 
tration of local f moments can be continuously varied 
between 1 (YbAl3 is a heavy fermion with the Kondo 
temperature of 600-700 K) and by replacing Yb with 
open f-shell with Lu with closed f-shell (LUAI3 is a nor- 
mal metal). Dilution by doping with non-magnetic ions 
should give rise to the disappearance of the long range 
order, and the associated HG at some critical doping. 
Moreover, in YbAl3 the low-temperature anomalies ob- 
served in susceptibility and specific heat suggest the pres- 
ence of a second energy scale of the order of 40 The 
observed anomalies are quickly suppressed by doping 13 
suggesting that this energy scale is related to the onset 
of Fermi liquid coherence i 12 i 13 Using optical pump-probe 
spectroscopy in the low perturbation regime we show that 
at high T and doping levels < x < 0.3 the relaxation 
rate decreases with decreasing T much like in other heavy 
electron systems studied so far.^ 5 -^ This behavior is con- 
sistent with the relaxation bottleneck due to the presence 
of the HG4 The HG does not change significantly in the 
doping range < x < 0.3, and the disappearance of the 
long range order appears near x ~ 0.4. 

At some characteristic doping dependent temperature 
T*(x), the slowing down of relaxation is truncated and 
the relaxation rate becomes constant below T* . This be- 
havior is explained by the competition between the relax- 
ation across the HG and the ballistic electron transport, 
when at low T the electron mean free path I becomes 
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larger than the optical penetration depth A. We derived 
the theoretical model that describes the functional form 
of relaxation for such a case and discuss the striking dop- 
ing dependence of the low temperature relaxation rate. 



II. EXPERIMENTAL 

The experiments were performed in the stan- 
dard pump-probe configuration using a mode-locked 
Ti:sapphirc laser as the source for both photoexcitation 
and probe laser pulse trains. The pump and probe pulses 
had a duration of < 20 fs at the center wavelength of 800 
nm. The experiments were performed in the very weak 
excitation regime with the excitation fluence of F w 0.4 
/iJ/cm 2 ^ From the literature values^ of the complex 
conductivity at 1.5 eV we estimated the corresponding 
absorbed energy densities to be w 90 mJ/cm 3 in YbAl 3 . 
Based on the known values of the Sommerfeld constant 
we estimated the corresponding increase in the electronic 
temperature after excitation to be less than 10 K over the 
entire T range implying that the dynamics being studied 
is near the thermal equilibrium. The Ybi-^Lu^Abj sam- 
ples (10 doping levels were studied) were grown by the 
"self-flux" method in excess Al.— 



III. RESULTS 

Figure 1 shows the dynamics of photoinduced reflectiv- 
ity change in YbAl3 and Ybo.95Luo.osAl3 in the T range 
between 4 K and 300 K. At 300 K, the rise time of about 
70 fs is followed by fast electronic relaxation which can 
be well described by a single exponential decay with the 
relaxation time r = 330 fs. Following this initial elec- 
tronic relaxation the system recovers to equilibrium on 
the timescale of several 100 ps, which is attributed to 
heat flow out of the excitation volume. Upon lowering 
T, changes in both the amplitude of the transient, A, 
and t are observed. Upon cooling t gradually increases 
in the same fashion as in previously studied heavy elec- 
tron compounds^^ reaching 3.2 ps at 50 K. There 
the anomalous increase in t is truncated and r remains 
constant below « 50 K. The amplitude of transient, A, 
also shows a pronounced T dependence similar to early 
reports^ At low temperatures A is being nearly constant 
up to ss 50 K and showing a pronounced suppression 
upon further T increase. Qualitatively similar results are 
obtained for all < x < 0.3, except that the temperature 
where the relaxation time divergence is truncated, T*, is 
doping dependent, as discussed below. In LuAl 3 , how- 
ever, the dynamics follows the behavior seen in normal 
metals, where the relaxation is only very weakly temper- 
ature dependent with the relaxation time slightly increas- 
ing upon increasing the temperature?^ Similarly, virtu- 
ally no T dependence of A is observed. 

The T dependence of r in YbAl 3 is shown in Figure 2. 
Above « 50 K it is consistent with the relaxation bottle- 
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FIG. 1: (color online) Photoinduced reflectivity traces in a) 
YbAls at selected temperatures and b) in Ybo.95Luo.05 AI3 
from 6-280 K (in 6K steps between 6-140K and 10K steps 
from 140-280 K). 

neck scenario due to the presence of the narrow HG, Eh g , 
in the density of statesi^ The relaxation process can be 
described by the RT model. Here following photoexcita- 
tion and the initial electron-electron and electron-phonon 
collisions which proceed on a sub-ps timescale, the sys- 
tem is characterized by excess densities of electron-hole 
pairs (EHPs) and high frequency phonons (HFP). When 
an EHP with an energy > Ehg recombines, a high fre- 
quency phonon (lo > Eh g ) is created. Since HFP can 
subsequently excite EHP, the recovery is essentially gov- 
erned by the decay of the HFP populations^ In the low 
excitation limit, as in this case, the T dependencies of 
both amplitude of the induced reflectivity change A, and 
t are governed by the T dependence of the number den- 
sity of thermally excited EHPs, nr- It was shown^i that 
n T cx A' 1 - 1, where A(T) = A(T)/A(T — > 0) and 

r(T)oc[5A(T) + 2n T (T)r 1 , (1) 

where S is a constant that depends only on the photoex- 
citation intensity. In a narrow band semiconductor tit 
depends on the shape of the DOS in the energy range 
e rj T around the chemical potential (here and further 
ks = 1). Generally, the number density of thermally 
excited EHPs across the gap (in this case indirect hy- 
bridization gap, Ehg) is given by 

nr*iT*exp(-E hg /2T) , (2) 

where p is of the order of 1 (0.5 for a BCS superconduc- 
tor), depending on the shape of the DOS near the gap 
edge. Neither the size of Eh g , the exact shape of the low 
energy DOS (anisotropy, possible impurity levels within 
the gap), or the T dependence of Eh g is well known in 
intermediate valence systems, therefore there is sonic am- 
biguity in determining Ehg- However, as the main T de- 
pendence in tit comes from the exponential term, a rough 
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estimate of the size of Eh g can be obtained. For T > 50 
K a good agreement of the data with this simple model 
can be obtained. Since r was found to be independent 
on F over the range of F studied^ we can assume that 
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Indeed Fig. 2 shows that there is a good 



agreement between r(T) and n T 1 (T) oc (A 1 — 1) 1 . Fit 
with n^ 1 cx T~ 5 exp(Eh g /2T) gives a values for the in- 
direct HG of Ei ig sa 16 meV. This value is lower than 
the values of the pseudogap obtained from the inelastic 
neutron scattering data (30 meV)^! as well as the value 
inferred from the optical conductivity data on YbAl3 (60 
meV).— The relationship between the observed pseudo- 
gaps in the spin (neutrons) and charge (optical) exci- 
tation spectra is unclear. On the other hand, the in- 
terpretation of the optical conductivity data is also not 
straightforward,-^ since the indirect transitions are for- 
bidden by the momentum conservation law. Further- 
more, the relaxation dynamics bottleneck is expected to 
be governed by the gap minimum, therefore taking the 
uncertainties in the interpretation of both experimental 
results the agreement is reasonable. Given the fact, that 
the mid infrared peak at 0.25 eV,— which is in the in- 
termediate valence systems commonly interpreted as the 
direct optical transition across the hybridization gap^ is 
in YbAl3 present way above room temperature^ it is 
expected that the same is true also for the indirect gap, 
as implied by the present study. 
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t(T > T*) can be well understood in terms of the RT 
scenario^ the results imply that the HG persists over 
wide range of doping. We should note that there is some 
variation in the gap extracted from the fit to the relax- 
ation times of individual doping levels (10 meV < Eh g < 
21 meV) - see the two fits in Figure 3. However, there is 
no systematic dependence of Ef lg (x). At x = 0.4, how- 
ever, a qualitative change in the relaxation dynamics is 
observed, with relaxation showing a very weak tempera- 
ture dependence, much like in the metallic LUAI3. This 
observation can be attributed to the loss of the long range 
order at x — 0.4, consistent with the inelastic neutron 
scattering data. 17 Based on the above we argue that the 
hybridization gap in Ybi-^Lu^Ala persists, nearly un- 
changed, up to x ~ 0.3, with the loss of the long range 
order to appear only near x ~ 0.4. 
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FIG. 3: (color online) The temperature dependence of t, ex- 
tracted using a single exponential decay fit, for Ybi-zLu^Als 
series. The dashed and dotted lines present fits with r oc 
T-°- 5 exp(E hg /2T) with E hg = 17 meV and E hg = 10 meV, 
respectively. 



FIG. 2: T dependence of r and the inverse density of ther- 
mally excited quasiparticles, n^ 1 , for YDAI3. n^ 1 was deter- 
mined from the T dependence of the normalized amplitude 
of the photoinduced reflectivity change, A (inset), via tit oc 
A~ L - 1. The data are fit with n T ~ T 1/2 exp(-E hg /2T), 
where Eh g = 16 meV is the size of the indirect HG. 

Figure 3 presents the T dependence of r, obtained by 
the single exponential fit, on a series of Ybi_ x Lua;Al3 
samples. Over a large range of doping, < x < 0.3, 
all t(T > T*) fall almost on the same curve. Since 



Below T*(x), however, r saturates at t*(x). Generally, 
the saturation of r at low T is expected within the RT 
model when tit{T) becomes smaller than 8A{T) - see 
Eq.([T]). However in this case the low temperature relax- 
ation time, r*, should depend on the excitation fluence. 
The fact that relaxation dynamics in this low perturba- 
tion experiment is found to be independent of excitation 
intensity effectively rules out this interpretation^ Alter- 
natively, we can assume (quite generally), that the relax- 
ation of the excited state proceeds via two independent 
channels. In the case of two competing relaxation pro- 
cesses the recovery dynamics will always be governed by 
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the fastest of the two. If one of the processes is slow and 
temperature independent, while the other shows strong 
temperature dependence, becoming slower and slower as 
the temperature is reduced, one would expect to ob- 
serve a crossover from the temperature dependent to a 
temperature independent relaxation at the temperature 
where the two timescales become comparable. Clearly, 
this scenario is consistent with the experimentally ob- 
served temperature dependence of the relaxation time. 
While at temperatures above T* the temperature depen- 
dence of both amplitude and relaxation time are in very 
good agreement with the relaxation bottleneck model - 
see Fig. (2), where the bottleneck in relaxation is a result 
of the presence of the hybridization gap near the Fermi 
level, the question arises which temperature independent 
relaxation process could lead to the observed crossover 
into T-independent relaxation below T*? 



IV. ANALYSIS AND DISCUSSION 

Measurements of the de Haas-van Alphen (dHvA) ef- 
fect in YbAl3 and LUAI3 indicate a large mean free 
path, I w 120 — 150 nni ) 19 i 20 several times larger than 
the optical penetration depth A (XybAh = 22 nm and 
XluAi 3 =17 nm)J£ From the T dependence of resistivity 
(inset to Fig. 5) it follows that / > A up to sa 50 K, 
well into experimentally accessible T range. Hence one 
can expect that in addition to the recombination pro- 
cess with characteristic time r a ballistic transport of 
hot electrons^ 1 - will take place on a timescale of X/vf, 
where vf is the Fermi velocity. From the dHvA dat a 19 i 20 
one obtains v F « 4 10 4 (1 10 6 ) m/s in YbAl 3 (LuAl 3 ) 
giving the characteristic timescale for the ballistic trans- 
port out of the probed volume of X/vf = 0.5 ps and 20 
fs, respectively. Since up is a weak function of tempera- 
ture, the competition between the strongly T dependent 
recombination across the HG and a T-independent bal- 
listic transport can account for the observed temperature 
dependence of the relaxation time in Ybi_ x Lu x Al3. 

To describe the relaxation process for such a case we 
use the Boltzmann kinetic equation with the collision in- 
tegral in the r-approximationj 2 ^ 



by the initial condition f(t — 0,x,e,0) = F(e,x), de- 
scribing the initial distribution of the hot quasiparticles 
after photoexcitation, has the solution: 

f(t, x, e, 9) = exp (-t/r(e, 0))F(e, x - v F cos (6,t)). (4) 

The spatial and time distribution of the electron- hole pair 
density is found as 

/>00 />7T 

n{x,t)= de f(t,x,e,9)sm{9)d9. (5) 
Jo Jo 

For comparison with the experimental data we assume 
that r is energy and angle independent, r(e, 9) — r and 
take the Gaussian form of the energy integrated excita- 
tion profile at t = 0, J °° F(e, x)de — exp (— x 2 /X 2 ) 

AR/R is proportional to the number of photoinduced 
carriers within A. Integrating Eq.([5]) with the Gaussian 
probe profile, exp (— x 2 /X 2 )^ and normalizing to 1 at 
t = one obtains 

AR/R cx n(t) = . -—e^^erftvpt/V^X). (6) 

V 2 VFt 

Eq.© has two regimes. On a short time scale, t < X/vp, 
it reduces to 
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n(t) = exp (-t/r - v 2 F t 2 /2X 2 ) 
while for t > X/vp we have 

A 



n(t) 



,-t/r 



V F t 



(7) 



(8) 



or 
2 



- + , F cos W - = --^ y . 




(3) 
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Here f(t,x,e,9) is the non-equilibrium correction to the 
equilibrium distribution function, which depends on time 
t, distance from the surface x, energy e relative to the 
Fermi energy, and the angle 9 between the velocity and 
the transport direction, x. The relaxation time r in gen- 
eral depends on the energy and the angle. Here we as- 
sume that the particle-hole symmetry is preserved, so 
that there is no electric field. Also, if the Fermi energy 
is large compared to the photon energy, one can neglect 
the dependence of the speed of hot electrons and holes on 
their relative energy, i.e. v w v F . Eq.Q, supplemented 



FIG. 4: (color online) The data recorded on YbAl3 at selected 
temperatures fit with single exponential decay (dashed blue 
line) and with Eq.© (solid red line), where X/vf = 1-33 ps. 
At low temperatures the fit is entirely governed by the X/vf 
term (with r > 5 ps), while above 50 K the relaxation is best 
fit with a single exponential decay. 

Figure 4 presents the photoinduced reflectivity traces 
in YbAl3 at a few selected temperatures. We fit the pho- 
toinduced reflectivity traces with Eq. © - solid lines and 
compare it to the single exponential decay (dashed). 
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At T < T* the model clearly describes the experimen- 
tal data much better than a single exponential decay. In 
fact, the fit is mainly governed by the X/vf ratio (best fit 
gives X/vp = 1-33 ± 0.05 ps in good agreement with the 
above estimate) while r can be anywhere between 5-100 
ps. The situation is reversed at T > T*, where the fit is 
entirely governed by r. Since the mean free path quickly 
decreases below A upon increasing T, e.g. in YbAl3 this 
should happen at w 50 K (see inset to Fig. 5), the ballis- 
tic transport becomes ineffective. Correspondingly, the 
fit to the 70 K data in Fig. 4 with fixed X/v F = 1.33 
clearly becomes inadequate. In fact, the data are much 
better fit with the single exponential decay (dashed blue 
line). Thus, T* can be thought of as the temperature 
where I ~ A; below T* the relaxation is dominated by 
the ballistic transport, while above this temperature re- 
laxation is governed by the recombination across the in- 
direct HG. In earlier studies^ such a saturation of r at 
low T was not observed. This can be understood, since 
at experimentally accessible temperatures of T > 10 K 
the mean free path in YbAl3 is by far the largest among 
the systems studied. Still, comparison of dynamics in 
YbCdCu 4 and YbAgCu 4 (Fig. 5 of Ref. [4]) shows an 
onset of saturation at w 20 K in YbAgCu4. This fea- 
ture is absent in YbCdCu4. Since the mean free path 
in YbCdCu4 is considerably smaller than in YbAgCu4, 25 
this observation is consistent with the proposed scenario. 
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FIG. 5: (color online) Doping dependence of t^ 1 for 
Ybi_a;Lu^Al3. Inset shows the T dependence of resistivity 
for several doping levels. 

The doping dependence of the low temperature relax- 
ation rate, t# (x), is plotted in Fig. 5. Initially t~ 1 
shows a weak decrease from t~ (a; = 0) « 0.2 ps -1 , fol- 
lowed by a rapid drop at x = 0.1. Further increase in 
x results in an increase in t^ 1 up to w 5 ps -1 . Within 
the above scenario, following the X/vp ratio, one would 
expect r~ 1 (LuAl 3 ) « 20 r7 1 (YbAl 3 ), as observed. The 



strong drop in r" 1 near x — 0.1 can be understood as 
being a result of the decrease of I below A, consistent 
with the resistivity data (see insert to Fig. 5). Indeed, 
in Ybo.9Luo.1Al3 the expected low T divergence of t due 
to the relaxation bottleneck^ is largely recovered (see 
Fig. 3). Thus, the qualitative doping dependence of r" 1 
in Ybi-zLuzAls can be understood as a competition be- 
tween decreasing X/vp ratio and doping induced disorder. 
Still, the underlying origin of the very sharp anomaly at 
x as 0.1 is unclear. Noteworthy, this doping level corre- 
sponds well with the doping level where low temperature 
coherence effects are strongly suppressedj 12 ' 13 



V. CONCLUSIONS 

We studied the carrier relaxation dynamics in 
Ybi-zLuzAbj, an intermediate valence heavy fermion 
system where doping was continuously tuned between a 
heavy fermion and normal metal. In YbAl3 there seems 
to be an overall agreement that its low temperature prop- 
erties are governed by a hybridization of conduction elec- 
trons with spatially extended wave functions and local- 
ized f orbitals, which can be described by the Anderson 
lattice model. Since photoexcited carrier relaxation was 
found to be very sensitive to the appearance of the gap in 
the charge excitation spectrum, the aim of this work was 
to utilize time-resolved optical spectroscopy to determine 
the doping level at which long range order, leading to the 
indirect hybridization gap in intermediate valence heavy 
electron systems, is suppressed. By performing system- 
atic studies of photoexcited carrier relaxation dynamics 
in Ybi-zLuzAls system, where doping was continuously 
varied between < x < 1, we show that the hybridiza- 
tion gap persists up to the critical doping x ~ 0.4. This 
is consistent with inelastic neutron scattering data. 17 In- 
deed, the percolation threshold for the cubic system is 
« 0.3i2£ This observation implies that the Kondo lattice 
is a very robust feature, and suggest the Anderson lattice 
model to be a minimum model for the description of this 
class of materials. 

In Ybi-aXuzAbj, below some doping dependent tem- 
perature T*(x), the relaxation time divergence due to 
the relaxation bottleneck is cut. This observation can be 
naturally explained by considering a parallel temperature 
independent relaxation channel. Since YbAl3 is known 
for its large electronic mean free path, we consider bal- 
listic electronic transport out of the probed volume as a 
possible competing process. We developed a theoretical 
model, describing the dynamics of excess quasiparticle 
density for such a case. Comparison with the experi- 
mental data showed that the functional form of relax- 
ation below T*(x), the temperature dependence of re- 
laxation time, as well as the doping dependence of the 
low temperature relaxation rate in Ybi-zLuz AI3 can be 
well accounted for by this model with the dominant bal- 
listic electron transport at low T. The strikingly sharp 
anomaly in t~ 1 {x) at x ~ 0.1 clearly requires further 



6 



experimental and theoretical work. ftskolleg and Center for Applied Photonics at the Uni. 

Konstanz, the Laboratory Directed Research and De- 
velopment program at Los Alamos National Laboratory, 
Acknowledgments and Center for Integrated Nanotechnologies at LANL. 

This work was supported by Sofja-Kovalevskaja Award 
from the Alexander von Humboldt Foundation, Zukun- 



1 A.C. Hewson, The Kondo Problem to Heavy Fermions, 
(Cambridge University Press, Cambridge, 1993). 

2 L. Degiorgi, Rev. Mod. Phys. 71, 687 (1999). 

3 PS. Riseborough, Adv. in Phys. 49, 257 (2000). 

4 J. Demsar, V.K. Thorsm0lle, J.L. Sarrao, and A.J. Taylor, 
Phys. Rev. Lett. 96, 037401 (2006) 

5 J. Demsar, J.L. Sarrao, and A.J. Taylor, J. Phys.-.Cond. 
Mat. 18 R281 (2006) 

6 J. Demsar, et al., Phys. Rev. Lett. 91, 027401 (2003). 

7 K. H. Ahn, et al., Phys. Rev. B 69, 045114 (2004). 

8 K.S. Burch et al., Phys. Rev. Lett. 100, 026409 (2008). 

9 A. Rothwarf, B.N. Taylor, Phys. Rev. Lett. 19, 27 (1967). 

J. Demsar, et al., Phys. Rev. Lett. 91, 267002 (2003). 

1 V.V. Kabanov, J. Demsar, D. Mihailovic, Phys. Rev. Lett. 
95, 147002 (2005). 

2 A. L. Cornelius, et al., Phys. Rev. Lett. 88, 117201 (2002). 

3 E. D. Bauer, et al., Phys. Rev. B 69, 125102 (2004). 

4 R.H.M. Groeneveld, R. Sprik, A. Lagendijk, Phys. Rev. B 
51 11433 (1995), and the references therein. 

5 Fluence was kept low to avoid sample heating. Over the 
fluence range F = 0.1 — 2 /ij/cm 2 no nonlinearities were 
observed, i.e. the photoinduced transient scales with F. 

6 S.J. Lee, et al., Phys. Rev. B 67, 075104 (2003). 



17 R. Osborn, et al., J. Appl. Phys. 85, 5344 (1999). 

18 H. Okamura, F. Michizawa, T. Nanba, T. Ebihari, J. Phys. 
Soc. Jap. 73, 2045 (2004); H. Okamura, et al., J. Mag. and 
Mag. Mat. 272-276, e51 (2004). 

19 T. Ebihara et al., J. Phys. Soc. Jap. 69, 895 (2000). 

20 I. Sakamoto et al., J. of Alloys and Compounds 323-324, 
623 (2001). 

21 see e.g. E. Knoesel, A. Hotzel, M. Wolf, Phys. Rev. B 57, 
12812 (1998). 

22 Such choice of the excitation profile is used to satisfy the 
condition dF/dx — at x = 0. This condition means 
that quasiparticles cannot escape from the sample through 
the surface. The results are independent of the particular 
choice of excitation profile. 

23 Analysis of the relaxation beyond the r-approximation 24 
provides about the same qualitatively picture (V. V. Ka- 
banov et al. in preparation). 

24 V.V. Kabanov, A.S. Alexandrov, Phys. Rev. B 78, 174514 
(2008). 

25 J.L. Sarrao et al., Phys. Rev. B 59, 6855 (1999). 

26 see e.g. CD. Lorenz, R.M. Ziff, J. of Physics A 31: 8147 
(1998). 



